The Balloon-borne Large Aperture Submillimeter Telescope 
(BLAST) 2005: Calibration and Targeted Sources 
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ABSTRACT 

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully 
during a 100-hour flight from northern Sweden in June 2005 (BLAST05). As part of the cali- 
bration and pointing procedures, several compact sources were mapped, including solar system. 
Galactic, and extragalactic targets, specifically Pallas, CRL 2688, LDN 1014, IRAS 20126+4104, 
IRAS 21078-H5211, IRAS 21307-1-5049, IRAS 22134-1-5834, IRAS 2301H-6126, K3-50, W 75N, 
and Mrk 231. One additional source, Arp 220, was observed and used as our primary calibrator. 
Details of the overall BLAST05 calibration procedure are discussed here. The BLAST obser- 
vations of each compact source are described, flux densities and spectral energy distributions 
are reported, and these are compared with previous measurements at other wavelengths. The 
250, 350, and 500 /^m BLAST data can provide useful constraints to the amplitude and slope of 
the submillimeter continuum, which in turn may be useful for the improved calibration of other 
submillimeter instruments. 
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1. Introduction 

The June 2005 flight of the Balloon-borne Large 
Aperture Submillimeter Telescope (BLAST) in- 
corporated a 2-m spherical primary mirror and 
large-format bolometer arrays operating at 250, 
350, and 500 ^m. A complete description of 



the B LAST instrument is given in iPascale et al. 



(|2008[ ). The BLAST wavelengths sample the peak 
of the spectral energy distribution (SED) for cool 
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dust (~ 10-40 K) and are in a regime which is 
difhcult or impossible to access from even the 
best ground-based sites. As a result, BLAST has 
the ability to conduct unique Galactic and ex- 
tragalactic submillimeter surveys with arcminute 
resolution and high sensitivity. BLAST'S primary 
scientific motivations are to study the angular and 
redshift distribution and evolution of high-redshift 
star-forming galaxies and to identify proto-stcUar 
cores and the earliest stages of star formation 
within Galactic molecular clouds. 

BLAST conducted a 100-hour flight, launch- 
ing from northern Sweden on 2005 June 12, and 
landing in northern Canada on 2005 June 16 
(BLAST05). One relatively shallow extragalac- 
tic f ield and several large ('^ 8 deg^) Galactic fields 
(see lChapin et al.ll2008l : lHargrave et "al]|2008[ ) were 
mapped, from which a large number of Galac- 
tic pre-stellar sources were extracted at very high 
signal-to-noise ratio. The results obtained from 
the BLAST flight from Antarctica in December 
2006 (BLAST06) wiU be discussed in future arti- 
cles. 

BLAST05 continuum observations were made 
with three focal plane arrays consisting of 139, 
88, and 43 detectors at 250, 350, and 500 /im, re- 
spectively. The optics and detector layout in the 
cryostat make simultaneous measurements possi- 
ble by three arrays, having a common field of 
view of about 14' x 7'. The detectors are silicon- 
nitride micromesh ( "spider-web" ) bolometric de- 
tectors cou2led__with_2£A (maximum optical effi- 
ciency. Griffin et al.ll2002n spaced feedhorn arrays 



Spitzer Space Telescope (jYoung et al.ll2004[ ). 

Together these observations represent a broad 
sampling of bright submillimeter sources that were 
available during the flight, given t he strict visi- 



(jTurner et al.ll2001[ ). The arrays are prototypes 
of those developed for the S PIRE instrument on 
Herschel (| Griffin et al.ll2004l) . 

In this paper we report the observations of com- 
pact sources that were targeted by BLAST largely 
for the purposes of fiux and pointing calibration. 
These sources include Pallas, CRL 2688, K3-50, 
W 75N, Mrk 231, and Arp 220, which have all 
been well studied in the submillimeter and are 
commonly used as primary or secondary stan- 
dards. We also conducted several observations of 
known bright proto-stars that were not otherwise 
included in our wider Galactic plane survey fields, 
specifically IRAS 20126-1-4104, IRAS 21078-f5211, 
IRAS 21307-1-5049, IRAS 22134-^5834, and IRAS 23011- 
Finally we observed the "starless" dense core 
LDN 1014 that was recently studied with the 



bility constraints of the telescope pPasca le et al. 



20081 ). The planets Mars and Uranus, often used 



as submillimeter calibrators, were not visible dur- 
ing this fiight. Although observations of Saturn 
were attempted, there were concerns about gain 
stability because it was relatively close to the di- 
rection of the Sun, and it may also have suffered 
from saturation effects in the detectors, and hence 
we abandoned attempts to use these particular ob- 
servations. After carefully assessing the system- 
atic effects in the relevant parts of the SEDs of 
the other sources, we chose to single out Arp 220 
as our primary flux calibrator. Using Arp 220, 
we determine that the BLAST05 calibration un- 
certainties are 12%, 10%, and 8% in the 250, 350, 
and 500 //m bands, respectively. 

This paper is divided into three main sections. 
Section [2] outlines the basic reduction steps and 
characterization of BLAST05 data. Section [3] 
discusses the performance of the warm optics in 
BLAST05 and the impact of the degraded per- 
formance on the resolution and sensitivity of the 
experiment. Section [4] describes in detail the ab- 
solute calibration derived from the primary flux- 
calibrator Arp 220, while Section [5] summarises 
the BLAST flux-densities and the simple fits of a 
modified black-body emission model to the SEDs 
of each of the targeted sources from the BLAST05 
flight. 

2. Data Reduction 

We now describe the main steps in the data re- 
duction and calibration process. Raw data from 
BLAST consist of a set of bolometer timestreams 
in voltage units, sampled at 100 Hz. The raw 
bolometer data are first cleaned for post-flight 
analysis. The data are de-spiked and then de- 
convolved to remove the electronics fi l ters f rom 
the timestreams; see IPantachon et al.l ( 20081) for 
specifics of the cleaning and deconvolution per- 
formed on the BLAST data. The cleaned data 
are combined with a post-fiight pointing solution 
p ilPascale et al.l 120081 ) to make maps at each fre- 
quency, taking advantage of the multiple detec- 
tors, as well as significant scan cross-linking, to 



minimize striping due to instrumental drifts. 

The bolometers in each array are corrected for 
relative gains, or flat-fielded, so that meaningful 
multi-bolometer maps can be generated. The flat- 
field corrections are determined using individual 
maps made for each bolometer from a single point- 
source calibrator, in this case, CRL 2688 (the 
brightest point-like source observed) . These scans 
were designed such that a fully-sampled map can 
be generated from each bolometer individually. 
For each bolometer, the total flux from the point 
source is integrated and a flat-fielding cocfhcient 
for each bolometer is then calculated as the ratio of 
the flux from that bolometer and the flux from an 
arbitrarily chosen reference bolometer. All subse- 
quent maps are generated from timestreams that 
apply the flat-fielding coefficient appropriately to 
each bolometer. 

To calculate the flux density from a point 
source, we adopted a matched-filtering technique 
similar to that used to extract point sources from 
several re cent extra-galactic submillimeter sur - 
veys (e.g.. ICoppin et alllJOol IScott et aDl2006h . 
A point spread function (PSF) is generated by 
stacking and averaging several point sources from 
various maps in telescope coordinates (azimuth 
and elevation relative to the telescope). We area- 
normalize the PSF, P, such that 



^P{x,y)6x6y 



(1) 



where the double sum is over all pixels in the PSF 
map, and 5x and Sy are the angular dimensions 
of a pixel such that the units of P are sr~^. In 
all BLAST maps, square pixels are used so that 
Sx = Sy. BLAST maps are calibrated in sur- 
face brightness units [Jy sr~^]. A map of a point 
source, M{x,y), can be modeled by the normal- 
ized PSF centered over the source, scaled by its 
flux density, S, or SP{x,y) ~ M{x,y). If each 
pixel in M(x, y) has an uncertai nty a(x,y), calcu- 



lated in the map-making process (jPantachon et al. 
20081 ). then we can write x^ for the model, assum- 
ing independant noise: 



x' = 



E 



iSP{x,y)-M{x,y)Y 



a'^{x,y) 



(2) 



hood flux density, 

_ T..^y{M{x,y)P{x,y)la\x,y)) 



(3) 



Minimizing this x^ results in the maximum likeli- 



This technique uses our knowledge of the beam 
as a model for the shape of un-resolved point- 
sources in BLAST maps. Fitting the amplitude 
of this template is optimal in a S/N sense, and 
particularly important for measurements of faint 
sources. This method produces smaller measure- 
ment uncertainties than simple aperture photom- 
etry, as pixels near the peak signal are weighted 
more heavily than pixels in the wings of the bright- 
ness distribution (aperture photometry weights all 
pixels equally). 

This technique gives a low statistical uncer- 
tainty both because we have intrinsically high 
signal-to-noise data and also because knowledge of 
the beam shape is used to weight the uncertain- 
ties; however, errors in the PSF will bias the result. 
To check for a bias we used simple aperture pho- 
tometry as a second measure of flux densities for 
CRL 2688. The difference between the matched 
fllter and aperture photometry flux densities of 
CRL 2688 is less than 3%. Since the matched filter 
provides significantly smaller statistical errors, it 
is used to extract flux densities from all point-like 
sources in this paper. 

The PSF for each individual bolometer varies 
slightly across the array. The stacked PSF takes 
this effect into account for the bulk of the map 
over which each bolometer contributes approxi- 
mately equally to the combined signal. However, 
the extreme edges of the map are only sampled by 
a fraction of the bolometers and therefore exhibit 
different effective PSFs. We extract flux densi- 
ties from maps made with single bolometers using 
the stacked PSF to assess the error introduced by 
adopting this incorrect template, and find at most 
a 10% bias. The bias is maximized at the very 
edges of the map, but affects flux densities mea- 
sured to within 14' (width of the array) of the 
edge of the map along the scan direction (approx- 
imately aligned with R.A.), and 7' (height of the 
array) in the transverse direction (approximately 
aligned with Dec). None of the measurements 
discussed in this paper are affected by this bias. 

The flat-fielding process was repeated for 
CRL 2688 observations made during the middle 



and end of the flight to check for stabiHty. Al- 
though no bias or trend is seen, the coefficients for 
each individual bolometer vary by 5% rms. This 
effect averages down when multiple bolometers 
are used, and is negligible in full bolometer-array 
maps as used in this paper. The 5% scatter in 
flatfield coefficients appears to be dominated by 
measurement uncertainties; therefore this value 
can be considered a very robust upper-limit on 
the systematic drift in the relative calibration, and 
is insignificant compared to the total calibration 
error budget given in Table [TJ 

2.1. Responsivity Variations 

Variations in bolometer loading (due to chang- 
ing sky emission, for example) and bolometer 
base-plate temperature cause changes in bolome- 
ter responsivities. To characterize and correct for 
these changes, a calibration lamp ("cal-lamp") is 
located in the optics box at the center of the re - 
flective Lyot stop (see §3 of IPascale et al.l 120081 ). 



allowing for the measurement of these changes. 
The cal-lamp is pulsed once every 10-15 minutes 
during flight. The cal-lamp is o f the same de- 



sign as the one used for SPIRE (jHargrave et al. 
20061 ). It is designed to provide a stable signal 



over the flight, so any change in cal-lamp ampli- 
tude as measured by a bolometer is directly pro- 
portional to a change in that bolometer's respon- 
sivity. The signal baseline is removed by fitting a 
line to a 550 ms segment of data before and after 
the cal-lamp pulse. A template cal-lamp profile 
from a raw bolometer timestream is fit to every 
cal-lamp pulse in every bolometer. The template 
is chosen from a typical bolometer and it has been 
verified that the template fit amplitude is within 
2% of a square-wave fit. The amplitude of the 
fit is interpolated over time and inverted to gen- 
erate a multiplicative scaling cal-lamp timestream 
for each bolometer. Fluctuations in the respon- 
sivity are less than 8% rms, dominated by diur- 
nal variations due to differences in atmospheric 
loading and thermal emission from the telescope's 
changing temperatures. On time scales of a typ- 
ical source map the fluctuations are less than 2% 
rms. A constant responsivity timestream for each 
bolometer is generated by applying the multiplica- 
tive scaling cal-lamp timestream. The large-scale 
fluctuations are suppressed and the final respon- 
sivity variations are less than 2% rms. 



3. BLAST05 Warm Optics Performance 

Dming the BLAST05 flight, the warm optics 
(primary and secondary mirrors) did not perform 
within the specifications. The beams as designed 
were expected to be close to diffraction limited 
and approximately Gaussian, exp(— 0^/2(Tg), with 
FWHM = 2(21n2)i/2crB = 32", 48", and 64", at 
250, 350, and 500 /im, respectively. In this case 
the full width at half power of the beam, FWHP, 
is equal to the FWHM, and the beam solid an- 
gle can be characterized as 51 = 7rQ;p/(41n2) with 
ttD = FWHP. The point-spread functions in all 
three BLAST bands were measured from mul- 
tiple observations of the proto-planetary nebula 
CRL 2688. The in-flight beam-shapes are shown in 
Figure [H Noting that the central structure in each 
is on a scale comparable to the diffraction limit, 
the beams are clearly far from ideal, distributing 
considerable power into an outer ring (hexagon) 
of diameter - 200". 

With the pre-flight predicted Noise Equivalent 
Flux Density (NEED) of 200 mJys^/^ we would 
expect l-cr surface brightness fluctuations at the 
nominal resolutions of ~NEFD/r2 or 11, 4.7, and 
2.6 MJy sr-isi/2 at 250, 350, and 500 /zm, re- 
spectively. The measured sensitivities were in fact 
8.8, 4.8, and 2.7MJysr"i s^/^^ respectively, which 
shows that BLAST did not suffer a reduced sen- 
sitivity to extended diffuse sources, i.e. detectors 
worked as planned and optical efficiency was as 
planned. 

There is, however, a marked loss in sensitiv- 
ity for point sources, because of the larger-than- 
expected beam solid angle. The FWHP of the 
degraded beam is ^ 200" (somewhat smaller at 
500 /im) , given the power in the ring. The reduc- 
tion in point-source sensitivity is < (200/32)^ ~ 40 
at 250 /tm. This is a pessimistic estimate because 
there is considerable higher-resolution structure in 
the beam, which one takes advantage of in the 
PSE fitting described above. Alternatively, one 
can examine the histogram of pixel brightnesses 
contributing to the PSF, and working from the 
brightest, find the beam area that accounts for half 
the power. We find au = 186", 189", and 189" 
for 250, 350, and 500 /xm, respectively, and so the 
sensitivity is reduced by factors of approximately 
34, 16, and 9. Despite this reduction, the tar- 
geted sources described in this paper were still de- 
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Fig. 1. — Point Spread Functions (PSFs), provided by observations of CRL 2688 (^ 15. 2p . for each of the 
three wavebands from BLAST05, generated by stacking several point source maps in telescope coordinates. 
The small circles represent the expected diffraction limited FWHM for each of the wavebands (Pascale et al. 
2008). The PSFs are significantly more extended than expected, reducing sensitivity to point sources and 
the ability to distinguish crowded sources amidst fluctuating cirrus. 



tected with high signal to no ise. Further da t a from 
BLAST05 are discussed in IChapin et all (|2008l ). 
Hargrave et al.l ([2008) , and other forthcoming pa- 



pers. 

Although it is not critical to the analysis here, 
we have investigated the origin of the beam degra- 
dation. One possibility is misalignment, shifting, 
and tilting of the secondary and primary with re- 
spect to the nominal optical axis. The launch 
was sufficiently violent to cause the inner-frame 
locking mechanism to fail; the inner-frame (which 
houses the primary, secondary, and camera) vi- 
olently hit both hardware elevation stops, which 
could have damaged the telescope, including the 
secondary support mechanism and the primary it- 
self. From visual inspection of the images (and 
power spectra) there is clearly information in the 
beams near the diffraction limit, such as the cen- 
tral dip at 350 ^m, which is not present at the 
other two wavelengths, and the width of the outer- 
ring. While these structures are suggestive a prob- 
lem in the telescope, numerous experiments with 
the optical design software ZEMAXlj, in which 
all the optical elements were shifted and/or tilted 
from their nominal position, we were unable to 
reproduce the distinctive features in the beams at 
all the wavelengths simultaneously. There appears 



^http : //www ■ zemax ■ com/ | 



to have been an optical problem beyond out-of- 
focus or out-of-alignment optics. The hexagonal 
shape of the beam does suggest damage to the 
carbon-fiber primary itself, which was constructed 
from six panels (segments). In addition to ef- 
fects of the violent launch, another contributing 
factor could have been the light rain the payload 
endured before launch. Carbon fiber becomes par- 
ticularly weak in a high humidity environment, 
possibly leading to delamination of the six panels 
either before or during the launch. Furthermore, 
there was the usual freeze-thaw cycle on ascent, 
where the ambient temperature dropped as low 
as — 65°C. Further analysis of the PSF (including 
Fourier Transform) was not performed as phase- 
information is required for proper analysis; as the 
mirror was fully destroyed upon landing there is 
little use in debugging the system further. 

4. Astronomical Flux-calibration 

The primary scientific goals of the BLAST ex- 
periment demanded a flux calibration-accuracy 
of better than 10% in all three BLAST pass- 
bands. Achieving this was complicated by the 
variable sky-visibility of BLAST due to the un- 
stable projected-latitude of the telescope gondola 
during the flight, and the restrictions on visibil- 
ity due to the Sun and Moon avoidance criteria, 
the orientation of sun-shields and other baffling. 



and the elev ation range (25-60 



inner- frame (Pascale et al. 200 



of the gondola's 
. Consequently 
BLAST had only limited access to the calibration 
sources commonly used at submillimeter and far 
infrared (FIR) wavelengths. 

Since the ecliptic-plane was not visible during 
the BLAST05 flight, no absolute flux-calibration 
of BLAST could be determined from observations 
of Uranus or Mars, for which model SEDs are 
known to have systematic uncertainties < 5% at 
submillimeter wavelengths ( Griffin fc Ortonlll993t 
I Wright 2007) . The pre-flight strategy for achiev- 
ing a 10% calibration accuracy, recognising the 
above visibility constraints, therefore forced us to 
identify other Galactic and extragalactic sources 
that could act as primary and secondary calibra- 
tors, with the following requirements: (i) avail- 
ability throughout the flight; (ii) already con- 
sidered, in some cases, as secondary calibrators 
for ground-based sub-mm telescopes and FIR 
satellites; (iii) well-constrained SEDs in the FIR 
to mm-wavelength regime, enabling accurate in- 
terpolation of the band-averaged flux-densities 
at BLAST wavelengths; (iv) bright (> 1 Jy 
at 500 /xm) and compact sources (with respect 
to the BLAST beam-size, < 20") that reside 
in regions with minimal spatial-structure in the 
Galactic foregrounds or backgrounds, allowing 
accurate subtraction of any extended emission. 
Given these criteria we scheduled regular obser- 
vations throughout the flight of bright embed- 
ded protostellar-sources and compact H II regions 
within the W58 and Cygnus-X molecular-cloud 
complexes (including K3-50, W 75N, and DR21), a 
post-AGB star with planetary nebula (CRL 2688), 
the asteroid 2 Pallas, and two ultraluminous infra- 
red galaxies (ULIRGs), Arp 220 and Mrk 231. 

In the following sub-section we justify our se- 
lection of the ULIRG Arp 220 as our primary cal- 
ibration source for the BLAST05 flight. 



10" 



4.1. 



Arp 220 - a primary flux-calibrator 
for BLAST05 



Arp 220 (At^ Il966l ) is a well-studied ULIRG 
at redshift z — 0.018 with a FIR luminosity of 
approximately 2 x 10^^ Lq. The sub-millimeter 
emission is believed to be due to dust heated 
primarily by active star-formation (at a rate of 
^ 2OOM0/yr) and to a lesser extent by an embed- 
ded Active Galactic Nucleus (AGN). This cmis- 
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Fig. 2. — Spectral Energy Distribution (SED) 
of Arp 220, the absolute flux calibrator for 
BLAST05. The best-flt model (heavy sohd- 
line) was constrained by the published data 
(blue triangles, discussed in Section [iTTjl . and ex- 
cluded BLAST05 measurements. The grey-lines 
shows the 68% confidence interval, estimated from 
100 Monte-Carlo simulations, about the best-fit 
model. Black diamonds indicate the model predic- 
tions for BLAST05 at 250, 350 and 500 ^m. The 
uncertainties associated with these predictions are 
11%, 10%, and 8% in each band, respectively. The 
best-fit parameters of the single-temperature dust 
model are given in Tabled 



sion is confined to a compact (< 2" diameter) 
region, as traced through interferometric imag- 
ing, by the distribution of the 1.3 mm dust con- 
tinn um, and the molecu lar gas in the C0(2-l)- 
line (jScoville et al.lll997^ that fuels the starburst 
activity. Arp 220 resides far from the Galactic 
plane (6 — -1-53°) and extended foreground Galac- 
tic contamination a s seen in the 100 /im map 
(jSchlegel et al.lll998[ ) does not affect the measured 
SED. Arp 220 is therefore a point-like source 
as observed by BLAST, as well as most ground- 
based and space-borne FIR to millimeter tele- 
scopes, making the construction of the FIR-to- 
millimeter SED insensitive to the aperture-sizes 
of the photometric measurements (Figure [2]) . 



The quoted measurement errors for the ground- 
based 350-1 100 /im photometry of Arp 220 are 
< 10%, while the absolute calibration uncertain- 
ties are typically 10-20%. In the case of shorter- 
wavelength FI R measurements from ISOPHOT 
at 60-200 /xm (iKlaas et al.ll200l[ ) and 120-200 ^m 
|Spinogl io et al.ll2002l ). the absolute calibration of 
ISOPHOT, which is estimated to be - 30%, also 
dominates the total photometric accuracy. More 
precise photometry, however, at 60 and 100 /zm is 
available from t he IRAS Revised B right Galaxy 
Sample (RBGS. ISanders et al.ll2003h . with statis- 
tical uncertainties that have been quoted at the 
^0.1% level. The difference between the absolute 
calibration of the RBGS and the original IRAS 
Bright Galaxy Samples ranges from 5-25%, and 
our analysis therefore assumes an additional 10% 
absolute calibration uncertainty on the IRAS flux 
densities. As we show below, combining all the 
above data, and accounting for systematic uncer- 
tainties, we have been able to constrain a sim- 
ple SED model of thermal dust emission that al- 
lows interpolated band-averaged flux-densities of 
Arp 220 to be estimated across the BLAST band- 
passes with an acceptable ~ 10% uncertainty (Fig- 
ure H]). 

In order to predict the BLAST band-averaged 
fluxes of Arp 220, a single-temperature modified- 
blackbody SED of the form S^ = A{v/vQy^B^{T) 
was first fit to existing submillimeter and FIR ob- 
servations described previously. In the above ex- 
pression, 5*1/ is the flux density, A is the ampli- 
tude, /3 is the emissivity-index of the radiating 
dust-grains, vo is fixed at c/250/im, and B„{T) 
is the Planck function for a blackbody with tem- 
perature T. 

The smooth SED is fit to the data using x^ min- 
imization. In addition to the measurement errors 
quoted in §4.1) we have assumed an additional cor- 
related error of 5% for all of the photometry, since 
the bulk of these instruments were calibrated using 
the same Uranus SE D, which is known to h ave an 
uncertainty of --5% (JGriffin fc Ortonlll993[) . Fur- 
thermore, we have been careful to account for cor- 
related errors using the full data covariance ma- 
trix, 

x2 = (S-S)C-i(S~Sf, (4) 

where S is the vector of measured flux densities, 
and S are the model predictions. The diagonal el- 
ements of C give the variances, af, calculated as 



the quadrature sums of the photometric and cal- 
ibration uncertainties described above. The off- 
diagonal elements of C are calculated as follows. 
For each pair of measurements taken with the 
same group and same instrument, i and j, calibra- 
tion uncertainties are assumed to be completely 
correlated such that Cij — (JiO-p where a'^ is the 
standard deviation of the calibration error com- 
ponent for the measurement. Monte Carlo sim- 
ulations are performed to characterize the range 
of models that are consistent with the data. One 
thousand mock data sets are generated by adding 
realizations of correlated Gaussian noise described 
by C to the measured flux densities S. Each one 
of these data sets is then re- fit using Equation [H 
The best-fit model (T = 41.7K and (3 = 1.3, for 
the observed SED, and its uncertainties, are shown 
in Figure [21 

The band-averaged flux densities for each 
BLAST wavelength are obtained for each SED 
model fit as 



Sb 



T(v)S{v)du, 



(5) 



where we use tildes to represent model quan- 
tities. Here the normalized filter transmission 
profiles T{h') are measured empiri cally using a 



Fourier transform spectrometer (jPascale et al. 
20081 ). From these 1000 sets of band-averaged 



fluxes, we obtain estimates of variances of the cal- 
ibrator brightnesses in each band, (cf), as well 
as the correlated error between the fluxes in each 
band, {uiCTj). 

Maximum-likelihood estimates for the flux 
densities of all BLAST05 sources are obtained 
from equation ([3]) using un-calibrated maps of 
Arp 220 in units of Volts. Multiple observations 
of CRL 2688 were made throughout the flight, 
allowing us to estimate systematic gain varia- 
tions not traced by the internal cal-lamp pulses 
at the 3% level. Since this gain variation is signif- 
icantly lower than the calibrator flux-uncertainty 
estimated above, we calculate single calibration 
coefficients Ci in each band for the entire flight as 
the ratio of the band-averaged fluxes to the raw 
fluxes from our maps. Note that this calibration 
depends explicitly on the assumption of a smooth 
thermal SED that neglects the effects of molecular 
emission lines. 

This problem has been examined in detail for 



SCUBA 850 /im continuum observations o f a va - 



Friesen et all (|2005[ ) 



riety of Galactic objects 
conclude that the line contamination is typically 
< 10% in hot Galactic molecular cores, and a sim- 
ilar conclusion was reached by Ijohnstone et alj 
(1,2003) in a study of protostellar sources in Orion, 
although in some exceptional cases the line- 
contamination of the 850 /im continuum fluxes in 
the most energetic sources reached levels > 50%. 
Similarly, line-contamination, due exclusively to 
CO(3-2) emission, has been measured to be <10% 
in 850 /Ltm continuum fluxes, for '^SO^^ of th e 
45 local galaxies observed by lYao et alJ (J2003l ). 
Unfortunately there is almost no information 
on the strengths of possible line-contamination 
in the BLAST pass-bands, including redshifted 
molecular-lines at <^ 300 /im, and thus we assume 
that < 10% line-contamination is present in the 
BLAST data. 

Including all of these sources of uncertainty, the 
final percentage errors (8-12%) in the absolute 
flux calibration, and a Pearson correlation matrix 
showing the relationship between errors in the dif- 
ferent bands, are given in Table [1] The calibra- 
tion coefficients are used to convert raw (Voltage) 
maps into calibrated (Jy) maps. All astronomical 
analysis is performed on the calibrated maps. 

5. BLAST05 observations of bright sources 

In this section we summarise the calibrated 
BLAST observations of all the bright targets pre- 
sented in this paper, paying attention to the level 
of agreement between their calibrated BLAST 
fluxes and the predictions from an interpolation 
of their measured SEDs. Although these bright- 
sources were selected to provide checks on the con- 
sistency of the calibration throughout the flight, 
and were not considered to be part of the sci- 
entific program, we derive some of the physical 
properties (temperatures, FIR luminosities, dust 
masses) that can be obtained from a simple model 
of thermal dust-emission that has been fit to the 
improved spectral coverage of their SEDs (^^ 15.1.1^ 

In addition to conducting large-area surveys 
of Galactic fields un dergoing active starformation 
(|Chapinet al.ll2008[ ). BLAST05 observed a total 
of 9 bright, individual point-like sources (Pallas, 
CRL 2688, IRAS 20126+4104, IRAS 21078+5211, 
IRAS 21307+5049, IRAS 22134+5834, IRAS 23011+6126 



Mrk 231, and Arp 220), complemented by obser- 
vations of more extended star-forming structures 
towards W75 N in Cyg X and K3-50 in W 58, as 
part of the overall calibration strategy. Although 
not used explicitly as a calibration source, we also 
present in this paper the BLAST05 observations 
of the compact low-mass protostar L1014-IRS to- 
wards the nearby Galactic dark-cloud LDN 1014. 

All BLAST05 calibration targets are in the 
northern hemipshere (23° < Dec < 61°; Tabled]). 
BLAST therefore provides valuable new photo- 
metric data at intermediate wavelengths to com- 
plement those obtained in the FIR and the (sub- 
)niillimeter regime. 

With the advantage of simultaneous observa- 
tions at 250, 350, and 500 /im, BLAST measures 
the thermal emission of these Galactic and extra- 
galactic objects with a high flux-calibration accu- 
racy in a wavelength-regime that traces the spec- 
tral curvature at wavelengths longer than the spec- 
tral peak, and shorter than the Rayleigh-Jeans 
limit for cold dust (< 50 K). In terms of modeling 
the SEDs, this short-subniillimeter wavelength- 
regime sampled by BLAST is critical for unlock- 
ing the degeneracies between the effect of dust- 
temperature, opacity and emissivity-index on the 
overall FIR-to-millimeter emission. The SPIRE 
instrument, using similar bolometer-arrays and 
identical filters to BLAST, following t he future 
launch of Herschel (JGriffin et al.l 120041 ) . will ex- 
tend this work with greater sensitivity and spatial 
resolution. 

It is not the purpose of this paper to provide a 
detailed account of physical properties and nature 
of each bright-source observed in BLAST05. Thus 
we only discuss briefly in the following sub-sections 
the collective results btained from these BLAST 
observations at 250, 350, and 500 /im. 

5.1. Spectral Energy Distributions 

Flux densities at 250, 350, and 500/imfor all 
the above point-sources are extracted using the 
matched-filter technique (eq. [3] outhned in §[2]), 
while for extended or blended-sources at the res- 
olution of BLAST05 (e.g. towards K3-50 and 
W 75N), we convolve the PSF with a simple model 
which is then fit to the data. In the latter case the 
model consists of multiple Gaussian sources, where 
we allow the amplitude, width, and position of the 



Gaussians to vary as parameters of the fit. 

Since the BLAST filters have large spectral 
widths, and colors sampled by the filters are a 
strong function of temperature (especially for T < 
25 K), a correction must be made to calculate 
monochromatic flux densities. Either a color cor- 
rection may be applied to the band-averaged flux 
densities, or effective filter wavelengths for each 
object may be quoted. We choose the former, and 
quote effective flux densities at precisely 250, 350, 
and 500 /im. Once the SED has been fit to the 
BLAST data by minimizing x^, to obtain the best- 
fit dust temperature and emissivity-index, this is 
used to calculate the correction 



S(y) = 5b^ 



Sv 



(6) 



where S-q is the band-averaged BLAST flux mea- 
surement, S{v) is the SED model flux density 
(evaluated at 250, 350, and 500 /im), and ^b is 
the SED model band-averaged flux. The BLAST 
color-corrected flux-densities for all bright sources 
presented in this paper are given in Table [21 

5.1.1. Single-temperature cold dust models 

Incorporating the BLAST measurements into 
the existing FIR-mm wavelength {^ 50-2000 fjjn) 
SEDs for each source, we fit only a single- 
temperature dust component while recognising 
that any model that utilizes shorter-wavelength 
mid-IR data must naturally include thermal emis- 
sion from hot dust (>100K). This decision is 
justified in the context of deriving the band- 
averaged BLAST fluxes, which require only a 
smooth and accurate fit to the SED over the re- 
stricted wavelength- range of the BLAST filters 
( Pascale et al.l |2008[) . Although alternative and 
more complex models (e.g. two-temperature com- 
ponents) can be fit to the same FIR-mm data, 
the difference in the quality (e.g. minimised x^) 
of the fit is not significant when compared to the 
single-temperature model. 

As previous studies of the FIR-mm SEDS of 
local galaxies and ULIRGs have shown (e.g. Eales 
et al. 2000, Lisenfeld, Isaak & Hills 2000), rely- 
ing on limited photometric data, with typically 3 
to 5 measurements in the FIR (from IRAS and 
ISO) and (sub-) millimeter data (at 850 /im and 
1.2 mm data), these SEDs can be well-represented 



by a single dust temperature in the range of 
25-50 K. Similar conclusions have been reached 
from the studies of dust in the Galactic ISM. 
Even with the addit ion of ground- b ased 350 and 



BenfordI 11999*: 'Hunter et al 



450 /i n i data (e.g. . . , _ 

2000t iDunne fc Baled I2OOII: i^Bcelen et all 



20061 



Coppin et al.l 12008^ single temperature fits work 



fairly well, and there is typically no need to fit 
multiple temperature components to SEDs unless 
one also brings in data at A -C 100 /im. More 
complicated fits would only be demanded with 
higher-fidelity spectral measurements covering a 
much wider range of wavelengths. 

Given this, we show the mid-IR to millimeter- 
wavelength SEDs of all the bright compact sources 
observed by BLAST05 and their best-fit single- 
temperature models in Figs. [21-El and present the 
best- fit model-parameters in Table [31 

Adopting a distance to each source (^ 15. 8p . we 
determine the FIR luminosity and dust mass from 
the SED fits, recognising that the adoption of two 
or more temperature components would modify 
the dust-mass estimate. Determination of bolo- 
metric fluxes and dust masses fo l lows t he proce- 
dures described in IChapin et al.l (|2008f ). Uncer- 
tainties in quantities derived from the SED fits 
in Table [3l are obtained from Monte Carlos simi- 
lar to those described in § 14.11 which now include 
BLAST data and their correlated calibration un- 
certainties. For each mock data set the quantity 
in question is derived from the model fit. These 
values are then placed in histograms in order to 
extract means and 68% confidence intervals. Fur- 
ther details on these Monte Carlo simulations are 



provided bv lChapin et al.l (|2008l ) 



5.2. CRL 2688 

CRL 2688 is a post-asymptotic giant-branch 
star with a proto-planetary nebula ("The Egg" 
iPrice fc Murdock 198 3) , and it is considered an 
archetypal object with which to study the evolu- 
tion of AGB red giants into bipolar planetary neb- 
ulae. The distance to CRL 2688 is not well-known, 
but is estimated to be approximately 1.2 5 kpc 
(ICrampton et al.lll975HCohen fc Kuhilll977l) . De- 



spite residing in the Galactic plane, CRL 2688 ex- 
hibits compact submillimeter emission and suffi- 
cient contrast relative to its surroundings to pro- 
vide a comm o nly used secondary calibrator (e.g. 
Sandelll IT99I [Jenness et al.ll2002l ). The submil- 
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Fig. 3. — CRL 2688 spectral energy distribu- 
tion, showing BLAST data and best-fitting mod- 
els. Symbols and lines are as those described in 
Fig. H ISO LWS data between 43-194 /xm are 



and Omont et al.l 



presented by ICox et al. 
tak en from Sandelll ( 199 



UM^ 



19961). Other data ar e 



Jenness et al.l (|2002( ). 



limeter emission arises from a dust shell with an 



extent of approximately 5" (| Jenness et al.l 120021) . 
making it point-Hke for BLAST. CRL 2688 was 
visited regularly throughout the flight and thus 
the observations were useful for tracking the PSF 
shape (Fig. 2) and relative antenna gain varia- 
tions, making it our primary flat-fielding calibra- 
tor. 

Fainter surface-brightness shells around CRL 2688 
have also been observed at 120-180 /zir iwith ISO at 
a radial distance of ~ 150" and 300" (|Speck et al. 
I2OOO), although this has recentl y been dispute d 



with Spitzer MIPS observations (|Do et al.l [20051) . 
In any case, these shells would be large compared 
to the BLAST beam, and of uniform surface- 
brightness, such that any faint emission would 
be removed from the measured fluxes through our 
baseline subtraction. No long-term submillimctcr 
variability has been dete cted in submillimete r ob- 
servations with SCUBA ( Jenness et al.l 120021 ) and 
hence archival data can be combined to perform 
our calibrations. 



Fig. 4. — Pallas spectral energy distribution. 
Symbols and lines are as those described in Fig. [2l 



The mid-IR to FIR SED has been accurately 
measured by IS O LWS between 43 and 194 ^m 



i Cox et al.lll996l). A t submillimeter wavelengths, 
Jenness et al.l (1200 21) report the integrated fluxes 
derived from SCUBA maps at 450 and 850/im, 
over a 40" aperture which exceed the earlier single- 
beam peak fluxes of Sandell (1994). Combining 
these data with the new BLAST measurements 
(Figure [3]) we flnd an acceptable fit to a single- 
temperature modified black-body model, deriving 
a dust-temperature of 210 K and emissivity-index 
of 0.4 ± 0.2. The SED of CRL 2688 in the FIR- 
submillimeter regime traced by BLAST is not as 
accurately constrained as that of Arp 220 as there 
are no FIR measurments with either IRAS or 
Spitzer MIPS due to saturation. We derive un- 
certainties in the BLAST fluxes of CRL 2688 that 
range from -- 10% at 500 ^m to --20% at 250 ^m. 

5.3. 2 Pallas 

2 Pallas is the th ird largest asteroid, with a di 



1996h . It has 



ameter of 530 km ([Mitchell et al 
a well-determined, eccentric orbit with a semi- 
major axis between 2.1 and 3.4 AU. At the time 
of the BLAST observations, 2 Pallas was at a 
distance of 2.2 AU. With an angular size < 1", 
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Pallas is point-like to BLAST. Errors in the 
derived flux-densities of Pallas at BLAST wave- 
lengths are 5-10%, and are due mostly to uncer- 
tainties in the size, shape, and albedo of the as- 
teroid. The current best estimates for the phys- 
ical properties of Pallas come from other submil- 
limeter and IR flux density measurements inte- 
gr ated into the "thermo-phy sical" model (TPM) 
of iMiiller fc Lagerrod ^QM)- The TPM has an 
uncertainty of around 5-10% for integrated fluxes 
in the ISO LWS band. The absolute calibration 
of the model has a supposed accuracy of ~ 10% in 
the longer- wavelength BLAST bands, comparing 
predictions with measurements from the JCMT 
(850 /im) and CSO (350 ^^m) (T. MiiUer, private 
communication). Although this uncertainty is 
similar to that for the SED of Arp 220, Arp 220 
was adopted as the absolute flux calibrator over 
Pallas since it is brighter, and its SED does not 
vary over the duration of the flight. Changes in 
apparent brightness due to orbital effects are neg- 
ligible on the timescale of t he individual BLAST 
observations (|Mtilleii 120051 ). and the variations 
due to the 7.8 hour rotational period of 2 Pal- 
las are less than 7% during the entire BLAST 
flight (T. MiiUer, private communication). Al- 
though the TPM contains a full shape description, 
and spin model (enabling SED predictions at any 
epoch), the uncertainty in the absolute calibration 
of the model flux is comparable to the predicted 
maximal brightness variations during the BLAST 
flight. For this reason, in Table 2 we have simply 
averaged together our 4 sets of Pallas observations, 
and we plot only the BLAST photometry points 
in Figure HI with no attempt to relate these to 
other data at significantly different epochs. The 
TPM predicts Pallas fluxes of 13.2, 6.8, and 3.4 
Jy which are consistant with the BLAST05 mea- 
surments, differing by 14, 8.6, and 7.3% at 250, 
350, and 500 /im respectively. The submillimeter 
spectral-index (3 in the SED of Pallas, as measured 
by BLAST, is 0.2 ± 0.3, which is the shallowest of 
all the sources described in this paper. 

5.4. L1014-IRS 

LDN 1014 (JLvndd Il962l ). which was classi- 
fled as a "starless core" due to the lack of a 
FIR IRAS detec tion and no signat ure of a con- 
tinuum outflow (j Young et al.l 120041 ) . is a nearby 
dark cloud that contains one of the lowest lumi- 




Wavelength [^m] 

Fig. 5. — L1014-IRS spectral energy distribution. 
Symbols and lines are as those described in Fig. [2] 
The fit is poor, partly because a single temper- 
ature modified blackbody may not be the best 
model, but also because of resolution issues with 
some of the photometry. We have not included the 
lower 450 and 850 /^m data points, which clearly 
disagree with the remaining photometry, proba- 
bly because of resolution effec ts. Non-BLAST 
data a re from lVisser et al.l (|2002f) and lYoung et al. 
(l2004l) . 



nosity proto-stellar systems, L1014-IRS. More re- 
cently, a bipolar molecular o utflow has been dis- 



covered (JBourke et al.l 120051) 



The distance to 
LDN 1014 is not flrmly determined , with esti- 
mates of ~200pc ( |Huard et all lJ OOJ). 400-900 



(iMorita et all l2006h and < 500pc (jShirlev et al. 
2007|). The FIR-millimeter SED of L1014-IRS 



£C 



derived from a physical model of a circumstellar 
disk heated by the central protostellar object, and 
constrained by recent Spitzer observations from 3- 
70 /xm predict a str ong spectral-peak at ~350 fim 
(|Younget al.ll2004h . Since L1014-IRS was unre- 
solved by BLAST, we were unable to observe the 
predicted density-profile. The BLAST fiuxes at 
350 and 500 /xm, however, are in good agreement 
with the model, but suggest that the SED peaks 
at shorter wavelengths, < 250 /im. 
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Fig. 6. — Left: 250 /iin image of the K3-50 region, with the brighter source K3-50A to the south-west 
and K3-50C to the north-east. Spectral energy distributions for K3-50A (middle panel) and K3-50C (right 
pan el). Symbols and lines are a s those described in Fig. [2l Other data are taken from IClegg et al.l ( 19761 ) 
and I Wvnn- Williams et al.l ([l977|). The resolution of IRAS makes it impractical to extract short wavelength 
flux densities for K3-50C. 



5.5. K3-50 

K3-50 is a group of compact H II regions within 
the star-formin g complex W 58, a t a distance of 
about 8.5kpc (jPeeters et al.ll2002f ). Despite the 
non-optimal PSF, BLAST clearly resolves at least 
2 sources separated by —2', which can be ident i- 
fied as K3-50A and K3-50C (JHoward et al.lll996l ). 
Figure [6] shows the 250 /xm image. Fluxes for 
K3-50A and K3-50C are obtained by simultane- 
ously fitting to the BLAST data a model of two 
2-dimensional Gausssian sources convolved with 
the BLAST PSFs, whose sizes, positions, and am- 
plitudes are parameters in the fit. An alternative 
model of two point-sources, fixed at the position s 
of K3-50A and K3-50C (from' Howard et al.lll996l ). 
produces an indistinguishable solution for the in- 
ferred fluxes. 



Thompson et al.l ( 20061 ) report an integrated 
450 iim SCUBA flux of 256 ± 79 Jy for K3-50A, 
which is consistent with our 500 /im photometry, 
despite issues with chopping and the restricted 
map-size for SCUBA making the comparison po- 
tentially complicated. Although there is some sug- 
gestion in the literature that these sources may be 
variablcl, we fit their individual FIR-millimeter 
SEDs (Figure [S]) with single-temperature modified 



black-body emission. K3-50C appears to have a 
dust temperature of only 9K, substantially colder 
than that in K3-50A which has dust radiating at 
temperatures of ^40 K (Table [3]), in good agree- 
ment with previous estimates b ased on FIR KAO 
data (JThronson fc Harpeiill979| ). This differece in 
dust-temperature may be due to the fact that K3- 
50C is embedded more deeply in the molecular 
cloud than K3-50A, or it may reflect the differ- 
ent evolutionary stages of the obj ects, with K3- 



50A g enerally agreed to be younger (jHoward et al. 



19961) 



5.6. W 75N 



The W 75N (|Westerhoutl Il958h and DR21 
field contains a set of young protostar/compact 
H II regions located in the Cygnus X molec- 
ular cloud complex. T hese sources are at a 
distance of about 3kpc (Campbell et al.l Il982t 



i 



Pipenbrink fc Wendkeiill988l ). although other es- 
timates suggest a c loser distance of around 1.7 kpc 
(| Jakob et al.l 120071 ). The region has been exten- 
sively studied over a wide wavelength range, with 
many detailed spectroscopic (including OH and 

other maser lines) as well as continuum obser- 

\ 1 i — — Zk 

vations (see e.g. iDavis et al.M2007[ ). Despite its 

extended emission, W 75N has been suggested as 



^http : //www, jach . hawali ■ edu/ JCMT/continuim/calibration/sens/ p%t¥A'f^lfc^l¥ifrfkiKt^^.''t?-faiil' 



or submillimeter observa- 
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Fig. 7. — Top-le ft: 250 /im BLAST i mage of the W 75N region with green contours showing the distribution 
of IR extinction (jKumar et al.ll2007l ). We divide the BLAST emission into 3 sources identified with: W 75N 
to the north; DR21(0H) the northern part of the merged pair; and DR21, the southernmost source. The 
1' pixel-size in this map hides much of the structure in the BLAST PSF. SEDs show data and best-fitting 
models for W 75N (top-right), DR21 (OH ) (lower- left) , and DR21 (lower- r ight). Symbols and Hue s are a s 
tho se described in Figjl The data from lClegget all (|l976[ ). iDavis et all (|2007l ). iHelou fc Walked (|l988l ). 
and I Harvey et al.l (|l977[ ) are not adopted in the fits becau s e of t he smaller aperture used to determine those 
flux density measurements. The data from iHarvey et al.l ( 1977 ) suggest a temperature of 65 K for W 75N 
and 40 K fo r DR21 (OH), and th ese are used as constraints in our SED fits. Other data are from lClegg et al.l 
(Il976[ ) and lColome et all (Il99,5[ ). 



13 



tions (|Sandel]||200"3h . The W 75N region contains 
at least 3 point-like objects which are resolved 
by B LAST: W 75N; DR21 (JDownes fc Rinehart 
19661 ). and DR21(0H). Flux densities for W 75N 
are obtained by the matched filter technique, as 
described in §[2] Since DR21(0H) and DR21 over- 
lap, their flux densities are derived from a model 
of point sources convolved with the BLAST PSF 
(as described in §5.1), where the positions and 
fluxes are the free parameters. These BLAST 
data suggest that DR21 has a higher value of /3 
and a lower dust-temperature than DR21 (OH), 
in agreemen t with what has b een found in other 
studies fe.g. IJakob et al.ll2007l l. 

In addition to the detection of the protostel- 
lar cores within the H II regions, and as can 
be seen in Figure [71 the BLAST data also show 
significant extended surface-brightness structure 
(- lOOOMJysr-i at 250//m). This is highly- 
correlated with the near -IR-extinction map s to- 
wards the same region ([Kumar et al.l 120071 ). in- 
dicating that it is real, rather than being map- 
making artefacts. . 

5.7. Mrk 231 



Mrk 231 ( MarkariarJ [19691 1. also known as 
IRAS 12540-1-5708, is the most luminous infrared 
galaxy within the local volume out to a distance 
of 175 Mpc. Mrk 231 hosts a Seyfert 1 nucleus, 
although the FIR emission is known to be domi- 
nated by star-fo rmation with in a region which is 
of order 1' (see [Lipari et al.l[2005L and references 
therein). It is therefore point-like for BLAST, 
given the measured PSF (Fig. 1). Mrk 231 
has a well measured SED and has been used 
as a template for photometric redshift est i mates 
of extragalactic sour ces (e.g. [Klaas et al.[ [2001 : 



Aretxaga et al.l [20051). There are u seful far-IR 



data from IRAS and ISO, although iKlaas et al. 



(|200l[ ) describe the existing /gO photometry as 
"partly distorted and uncertain" . [BenfordI (|l999l ) 
reports a measurement at 350 /im of 1.05 ±0.25 Jy, 
which is consistent with the BLAST 350 /im pho- 
tometry. There is also 450 /xm photometry pub- 
lished, although using a much smaller aperture 



(JRigopoulou et all [l996[) . BLAST flux densities 



thus help constrain the SED in the poorly-sampled 
submillimeter regime (see Figure [HJ. However, 
given the relative faintness of this source, the im- 
provement is modest compared with other sources 




Wavelength [^m] 



Fig. 8. — Mrk 231 spectral energy distribution 
showing data and best-fitting model. Symbols and 
lines are as those described in Fig.[2l Mrk 231 was 
not detected by BLA S T at 5 iim. Other d ata ar e 
taken from [Benford[ ([l999l ). [Klaas et aL|j200l[ ). 



Lisenfeld et all ([2000D. [Stickel et al.[ ([2004D . and 
Helou fc Walked ([l988h . 



reported here. We find slightly lowe r values for 



tempe rature and /3 compared with [Klaas et al. 



([200 1[ ). but they also claim that the data are bet- 
ter fit with multiple temperature components. 

5.8. Compact Protostellar Sources 

The remaining Galactic targets, described be- 
low, are all compact protostellar sources, iden- 
tified via colour-criteria in the FIR, that satisfy 
the BLAST selection requirements. We present 
here only their fiuxes, dust masses and luminosi- 
ties (Tables 2 & 3) derived from the best-fit models 
describing their SEDs (Fig [^ . 

5.8.1. IRAS 20126+4104 

IRAS 20126-f 4104 is a high-mass compact pro- 
tostar, located towards the Cygnus X molecu- 
lar cloud complex, although its exact distance 
i s uncertain. We ad opt a distance of 1.7 kpc 
( Shepherd et al.[ [20001 ) . It is probably the best 
studied example of a massive proto-stellar system 
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Fig. 9.— 
Table [I]). 
a) IRAS 201 2 6+41 04 — Additional data are taken from iHelou fc Walkeij (I1988D. iHunter et al 



Spectral energy distributions of compact protostellar sources, including the BLAST fluxes (see 
The data symbols and fitted curves have the s ame meaning as those d escribed in Fig 



Beuther et al.l ([2002), and IWilliams et al 



(200i). The 350 ^m data point from 'H unter et all (|2o6oh is not 
used in the fit, d ue to the s maller aperture. - ( b) IRAS 21078+5211 — Additional data are taken from 



i). 



Helou fc Walked (|l988f ) and iJenness et al.l (|1995| ). Three data points are not used in the fit due to the 
smal l effective apert ure sizes in those me asureme nts. - (c) IRAS 21307+5049 — Additional data are taken 



from lMolinari et al.l (jJOOOi) and lHelou fc Walkeii (|l988l ). The 850 /im data point from iMolinari et al.l (|2000l ) 
is not used in the fit due to its s maller aperture compare d to the BLAST PSF. - (d) IRAS 22134+ 5834 



Additional data are taken from iHelou fc Walker 
(e) IRAS 23011+6126 



Williams et al. (20041 



and IChini et al.l (|200lh . The 850 /im data point from IChini et all (|200ll ) 



1988f) . IChini et al.l (|200l|), iBeuther et a T ('2002'), and 

Additional data are taken from lHelou fc Walker (19881 ) 

is not used in the fit due to the 



smaller aperture. 
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assoc iated with a disk and outflow (JEdris et al. 
20051 ). BLAST prov ides bet ter constrain ts on 



the SED (jWilking et a l. 1989; Hunter et al.l l20001. 
particularly the submillimeter slope, as shown in 
Figure [51 

5.8.2. IRAS 21078+5211 

IRAS 21078+5211 is a compact H II re- 
gion with a water maser, l ocated at a distance 
of approximately 1.65 k pc (JMolinari et al.l Il996t 



Wouterloot et all Il993h . 



The best-fit single- 
temperature model with T~34K and_^~l.l, is 



consis tent with a previous study by iMcCutcheon et al 
()l99lh . 



5.8.3. IRAS 21307+5049 



IRAS 21307-^5049 (also Mol 136. lMolinari et al. 

19961 ) is a high mass protostellar candidate. At a 
distance of 3.6 kpc, it has an angular size of about 



4" in the millimeter ( Molinari et ah 2002^. and 



is the refore point-like to BLAST. iFontani et al. 
(|2004l ) present 850^m SCUBA continuum observa- 
tions, detecting an extended halo of 40" diameter 
and a compact unresolved core. The extended 
halo is not detected by BLAST. We only fit a 
single-temperature modified blackbody curve to 
the cold core component. We find a temperature 
of 32 K and luminosit y of Lfir, — 4,000 L(d, in very 



good agreement with iMofinari et al.1 (|2000[ ) 



5.8.1 IRAS 22134+5834 

IRAS 22134-1-5834 is a medium-to-high mass 
protostar located within the extensive H II re- 
gion S134, in the constellation Ce pheus. It is at 



a dist ance of approximately 900 pc (jDobashi et al. 



19941 ). and has a total luminosity of ~ 1600 Lq. 



5.8.5. IRAS 23011+6126 

IRAS 2301H-6126 is a you ng protostar located 
at a distance of '^ 730 pc ([Crawford fc Barned 



19701). IRA S 2301 1H-6126 is point-like to BLAST. 
Chini et al.l ([20011 ) discusses the current best SED 
(similar to IRAS 22134-1-5834 above). The peak of 
the SED is at about 150 ^m, as shown in Figure[9l 
It appears to have T ~ 28 K, with a shallower slope 
(/3 = 0.6 ±0.2) than most of the other protostellar 
sources that BLAST observed. 



6. Conclusions 

For the BLAST 2005 flight, the out-of-focus 
PSFs significantly reduced the point-source sensi- 
tivity, and yet because this is such a difficult wave- 
band to observe from the ground, the BLAST05 
data still produced useful results. The submillime- 
ter fiux densities reported here agree well with 
other measurements, at least when the aperture 
sizes are consistent. However, BLAST05 uncer- 
tainties are typically much smaller than those of 
other measurements in the 250-500 /im regime. 

A particular advantage of the BLAST data is 
the ability to constrain the emissivity index /3, due 
primarily to the 3 separate submillimeter chan- 
nels, observed simultaneously with the same tele- 
scope, and having a consistent calibration proce- 
dure performed across all 3 bands. 

Arp 220 was adopted as the primary calibrator 
for the BLAST05 data. A useful cross-check is to 
note that the derived fiux densities for Pallas were 
all within la of the value s derived from the TPM 
(|Mtiller fc Lagerrodlioiol . 

By applying this careful calibration procedure, 
BLAST has been able to improve estimates of dust 
temperatures and masses for several relatively 
bright sources, and also to provide a database of 
submillimeter measurements which may be useful 
for calibrating future instruments. 

The BLAST collaboration acknowledges the 
support of NASA through grant numbers NAG5- 
12785, NAG5-13301, and NNG0-6GI11G, the 
Canadian Space Agency (CSA), Canada's Nat- 
ural Sciences and Engineering Research Council 
(NSERC), and the UK Particle Physics & As- 
tronomy Research Council (PPARC). We would 
also like to thank the Columbia Scientific Bal- 
loon Facility (CSBF) staff for their outstanding 
work, as well as T.G. Miiller for valuable discus- 
sions on Pallas. LO acknowledges partial support 
by the Puerto Rico Space Grant Consortium and 
by the Fondo Istitucional para la Investigacion of 
the University of Puerto Rico. CBN acknowledges 
support from the Canadian Institute for Advanced 
Research. This research has been enabled by the 
use of WestGrid computing resources. This re- 
search also made use of the SIMBAD database, 
operated at the Centre de Donees astronomiques 
de Strasbourg (CDS), Strasbourg, France, and 



16 



the NASA/IPAC Extragalactic Database (NED), 
which is operated by the Jet Propulsion Labo- 
ratory, Cahfornia Institute of Technology, under 
contract with the National Aeronautics and Space 
Administration. 

REFERENCES 

Aretxaga, I., Hughes, D. H., & Dunlop, J. S. 2005, 
MNRAS, 358, 1240 

Arp, H. 1966, ApJS, 14, 1 

Beelen, A., Cox, P., Bcnford, D. J., Dowell, C. D., 
Kovacs, A., Bertoldi, F., Omont, A., & Carilli, 
C. L. 2006, ApJ, 642, 694 

Benford, D. J. 1999, PhD thesis, Cahfornia Insti- 
tute of Technology 

Beuther, H., Schilke, P., Menten, K. M., Motte, F., 
Sridharan, T. K., & Wyrowski, F. 2002, ApJ, 
566, 945 

Bourke, T. L., Crapsi, A., Myers, P. C, Evans, II, 
N. J., Wilner, D. J., Huard, T. L., J0rgensen, 
J. K., & Young, C. H. 2005, ApJ, 633, L129 

Campbell, M. F., Niles, D., Nawfel, R., Hawry- 
lycz, M., Hoffmann, W. F., & Thronson, Jr., 
H. A. 1982, ApJ, 261, 550 

Chapin, E. et al. 2008, ApJ, in press 

Chini, R., Ward-Thompson, D., Kirk, J. M., Niel- 
bock, M., Reipurth, B., & Sievers, A. 2001, 
A&A, 369, 155 

Clegg, P. E., Rowan- Robinson, M., & Ade, 
P. A. R. 1976, AJ, 81, 399 

Cohen, M. & Kuhi, L. V. 1977, ApJ, 213, 79 

Colome, C, Harvey, P. M., Lester, D. F., Camp- 
beU, M. F., & Butner, H. M. 1995, ApJ, 447, 
236 

Coppin, K., Chapin, E. L., Mortier, A. M. J., 
Scott, S. E., Borys, C, Dunlop, J. S., Halpern, 
M., Hughes, D. H., Pope, A., Scott, D., Ser- 
jeant, S., Wagg, J., Alexander, D. M., Al- 
maini, O., Aretxaga, I., Babbedge, T., Best, 
P. N., Blain, A., Chapman, S., Clements, D. L., 
Crawford, M., Dunne, L., Eales, S. A., Edge, 
A. C, Farrah, D., Gaztaiiaga, E., Gear, W. K., 



Granato, G. L., Grevc, T. R., Fox, M., Ivison, 
R. J., Jarvis, M. J., Jcnness, T., Lacey, C, Lep- 
age, K., Mann, R. G., Marsden, G., Martinez- 
Sansigre, A., Oliver, S., Page, M. J., Peacock, 
J. A., Pearson, C. P., Percival, W. J., Priddey, 
R. S., Rawlings, S., Rowan-Robinson, M., Sav- 
age, R. S., Seigar, M., Sekiguchi, K., Silva, L., 
Simpson, C, Small, I., Stevens, J. A., Takagi, 
T., Vaccari, M., van Kampen, E., & Willott, 
C. J. 2006, MNRAS, 372, 1621 

Coppin, K., Halpern, M., Scott, D., Borys, C, 
Dunlop, J., Dunne, L., Ivison, R., Wagg, J., 
Aretxaga, I., Battistelli, E., Benson, A., Blain, 
A., Chapman, S., Clements, D., Dye, S., Far- 
rah, D., Hughes, D., Jenness, T., van Kampen, 
E., Lacey, C, Mortier, A., Pope, A., Priddey, 
R., Serjeant, S., Small, I., Stevens, J., & Vac- 
cari, M. 2008, MNRAS, 123 

Cox, P., Gonzalez- Alfonso, E., Barlow, M. J., Liu, 
X.-W., Lim, T., Swinyard, B. M., Cernicharo, 
J., Omont, A., Caux, E., Gry, C, Griffin, M. J., 
Baluteau, J.-P., Clegg, P. E., Sidher, S., Pe- 
quignot, D., Nguyen-Q-Rieu, King, K. J., Ade, 
P. A. R., Towlson, W. A., Emery, R. J., Furniss, 
I., Joubert, M., Skinner, C. J., Cohen, M., Ar- 
mand, C, Burgdorf, M., Eward, D., di Giorgio, 
A., Molinari, S., Texier, D., Trams, N., linger, 
S. J., Glencross, W. M., Lorenzetti, D., Nisini, 
B., Orfei, R., Saraceno, P., & Serra, G. 1996, 
A&A, 315, L265 

Crampton, D., Cowley, A. P., & Humphreys, 
R. M. 1975, ApJ, 198, L135 

Crawford, D. L. & Barnes, J. V. 1970, AJ, 75, 952 

Davis, C. J., Kumar, M. S. N., SandeU, G., Froe- 
brich, D., Smith, M. D., & Curric, M. J. 2007, 
MNRAS, 374, 29 

Do, T., Morris, M., Sahai, R., & Stapelfeldt, K. R. 
2005, in Bulletin of the American Astronomical 
Society, 1161 

Dobashi, K., Bernard, J.-P., Yonckura, Y., & 
Fukui, Y. 1994, ApJS, 95, 419 

Downes, D. & Rinehart, R. 1966, ApJ, 144, 937 

Dunne, L. & Eales, S. A. 2001, MNRAS, 327, 697 

Edris, K. A., Fuller, G. A., Cohen, R. J., & Etoka, 
S. 2005, A&A, 434, 213 



17 



Fontani, F., Cesaroni, R., Testi, L., Molinari, S., 
Zhang, Q., Brand, J., & Walmsley, C. M. 2004, 
A&A, 424, 179 

Friesen, R. K., Johnstone, D., Naylor, D. A., & 
Davis, G. R. 2005, MNRAS, 361, 460 

Griffin, M. J., Bock, J. J., & Gear, W. K. 2002, 
Apphed Optics, 41, 6543 

Griffin, M. J. & Orton, G. S. 1993, Icarus, 105, 
537 

Griffin, M. J., Swinyard, B. M., & Vigroux, L. 
2004, in Optical, Infrared, and MiUimeter Space 
Telescopes. Edited by Mather, John C. Pro- 
ceedings of the SPIE, Volume 5487, ed. J. C. 
Mather, 413-424 

Hargrave, P., Waskett, T., Lini, T., & Swinyard, 
B. 2006, in Millimeter and Submillimeter De- 
tectors and Instrumentation for Astronomy III. 
Edited by Zmuidzinas, Jonas; Holland, Wayne 
S.; Withington, Stafford; Duncan, William D.. 
Proceedings of the SPIE, Volume 6275, pp. 
627514 (2006). 

Hargrave, P. et al. 2008, ApJ, submitted 

Harvey, P. M., Campbell, M. F., & Hoffmann, 
W. F. 1977, ApJ, 211, 786 

Helou, G. & Walker, D. W., eds. 1988, Infrared 
astronomical satellite (IRAS) catalogs and at- 
lases 

Howard, E. M., Pipher, J. L., Forrest, W. J., & de 
Free, C. G. 1996, ApJ, 460, 744 

Huard, T. L., Myers, P. C., Murphy, D. C., Crews, 
L. J., Lada, C. J., Bourke, T. L., Crapsi, A., 
Evans, II, N. J., McCarthy, Jr., D. W., & 
Kulesa, C. 2006, ApJ, 640, 391 

Hunter, T. R., Churchwell, E., Watson, C, Cox, 
P., Benford, D. J., & Roelfsema, P. R. 2000, 
AJ, 119, 2711 

Jakob, H., Kramer, C, Simon, R., Schneider, N., 
Osscnkopf, v., Bontemps, S., Graf, U. U., & 
Stutzki, J. 2007, A&A, 461, 999 

Jenness, T., Scott, P. F., & Padman, R. 1995, MN- 
RAS, 276, 1024 



Jenness, T., Stevens, J. A., Archibald, E. N., 
Economou, F., Jessop, N. E., & Robson, E. I. 
2002, MNRAS, 336, 14 

Johnstone, D., Boonman, A. M. S., & van 
Dishoeck, E. F. 2003, A&A, 412, 157 

Klaas, U., Haas, M., Miiller, S. A. H., Chini, R., 
Schulz, B., Coulson, I., Hippelein, H., Wilke, 
K., Albrecht, M., & Lemke, D. 2001, A&A, 379, 
823 

Kumar, M. S. N., Davis, C. J., Grave, J. M. C, 
Ferreira, B., & Froebrich, D. 2007, MNRAS, 
374, 54 

Lipari, S., Terlevich, R., Zheng, W., Garcia- 
Lorenzo, B., Sanchez, S. F., & Bergmann, M. 

2005, MNRAS, 360, 416 

Lisenfeld, U., Isaak, K. G., & Hills, R. 2000, MN- 
RAS, 312, 433 

Lynds, B. T. 1962, ApJS, 7, 1 

Markarian, B. E. 1969, Astrofizika, 5, 286 

McCutcheon, W. H., Sato, T., Dewdncy, P. E., & 
Purton, C. R. 1991, AJ, 101, 1435 

MitcheU, D. L., Ostro, S. J., Hudson, R. S., 
Rosema, K. D., Campbell, D. B., Velez, R., 
Chandler, J. F., Shapiro, I. I., Giorgini, J. D., 
& Yeomans, D. K. 1996, Icarus, 124, 113 

Molinari, S., Brand, J., Cesaroni, R., & Palla, F. 
1996, A&A, 308, 573 

— . 2000, A&A, 355, 617 

Molinari, S., Testi, L., Rodriguez, L. F., & Zhang, 
Q. 2002, ApJ, 570, 758 

Morita, A., Watanabe, M., Sugitani, K., Itoh, 
Y., Uehara, M., Nagashima, C, Ebizuka, N., 
Hasegawa, T., Kinugasa, K., & Tamura, M. 

2006, PASJ, 58, L41 

Miiller, T. G. 2005, Private Communication 

Miiller, T. G. & Lagerros, J. S. V. 2002, A&A, 
381, 324 

Omont, A., Moseley, S. H., Cox, P., Glaccum, W., 
Casey, S., Forveille, T., Chan, K.-W., Szczerba, 
R., Loewenstein, R. F., Harvey, P. M., & Kwok, 
S. 1995, ApJ, 454, 819 



18 



Pantachon, G. et al. 2008, ApJ, in press 

Pascale, E. et al. 2008, ApJ, in press 

Peeters, E., Martin- Hernandez, N. L., Daniour, 
F., Cox, P., Roelfsema, P. R., Baluteau, J. -P., 
Tielens, A. G. G. M., Churchwell, E., Kessler, 
M. F., Mathis, J. S., Morisset, C., & Schaerer, 
D. 2002, A&A, 381, 571 

Pipenbrink, A. & Wendker, H. J. 1988, A&A, 191, 
313 

Price, S. D. & Murdock, T. L. 1983, AFGL-TR- 
0208 Environemental Research papers, 161, 1 

Rigopoulou, D., Lawrence, A., & Rowan- 
Robinson, M. 1996, MNRAS, 278, 1049 

Sandell, G. 1994, MNRAS, 271, 75 

Sandell, G. 2003, in ESA SP-481: The Calibra- 
tion Legacy of the ISO Mission, ed. L. Metcalfe, 
A. Salama, S. B. Peschke, & M. F. Kessler, 439- 
442 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, 
Surace, J. A., & Soifer, B. T. 2003, AJ, 126, 
1607 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 
1998, ApJ, 500, 525 

Scott, D., Chapin, E., Aretxaga, L, Austerniann, 
J., Coppin, K., Crowe, M., Frey, L., Gibb, A., 
Halpern, M., Hughes, D., Kang, Y., Kim, S., 
Lowenthal, J., Perera, T., Pope, A., Scott, K., 
Wilson, G., & Yun, M. 2006, in American As- 
tronomical Society Meeting Abstracts, 125.04 

ScoviUe, N. Z., Yun, M. S., & Bryant, P. M. 1997, 
ApJ, 484, 702 

Shepherd, D. S., Yu, K. C, Bally, J., & Testi, L. 
2000, ApJ, 535, 833 

Shirley, Y. L., Claussen, M. J., Bourke, T. L., 
Young, C. H., & Blake, G. A. 2007, ApJ, 667, 
329 

Speck, A. K., Meixner, M., & Knapp, G. R. 2000, 
ApJ, 545, L145 

Spinogho, L., Andreani, P., & Malkan, M. A. 2002, 
ApJ, 572, 105 



Stickel, M., Lemke, D., Klaas, U., Krause, O., & 
Egner, S. 2004, A&A, 422, 39 

Thompson, M. A., Hatchell, J., Walsh, A. J., Mac- 
Donald, G. H., & Millar, T. J. 2006, A&A, 453, 
1003 

Thronson, Jr., H. A. & Harper, D. A. 1979, ApJ, 
230, 133 

Turner, A. D., Bock, J. J., Beeman, J. W., Glenn, 
J., Hargrave, P. C, Hristov, V. V., Nguyen, 
H. T., Rahman, F., Sethuraman, S., & L., 
W. A. 2001, Appl. Opt., 40, 4921 

Visser, A. E., Richer, J. S., & Chandler, C. J. 

2002, AJ, 124, 2756 

Wester hout, G. 1958, Bull. Astron. Inst. Nether- 
lands, 14, 215 

Wilking, B. A., BlackweU, J. H., Mundy, L. G., & 
Howe, J. E. 1989, ApJ, 345, 257 

WiUiams, S. J., Fuller, G. A., & Sridharan, T. K. 
2004, A&A, 417, 115 

Wouterloot, J. G. A., Brand, J., & Fiegle, K. 1993, 
A&AS, 98, 589 

Wright, E. L. 2007, arXiv:astro-ph/0703640 

Wynn- Williams, C. G., Matthews, K., Werner, 
M. W., Becklin, E. E., & Neugebauer, G. 1977, 
MNRAS, 179, 255 

Yao, L., Seaquist, E. R., Kuno, N., & Dunne, L. 

2003, ApJ, 588, 771 

Young, C. H., J0rgensen, J. K., Shirley, Y. L., 
Kauffmann, J., Huard, T., Lai, S.-P., Lee, 
C. W., Crapsi, A., Bourke, T. L., DuUemond, 
C. P., Brooke, T. Y., Porras, A., Spiesman, 
W., Allen, L. E., Blake, G. A., Evans, H, 
N. J., Harvey, P. M., Koerner, D. W., Mundy, 
L. G., Myers, P. C, Padgett, D. L., Sargent, 
A. L, Stapelfeldt, K. R., van Dishoeck, E. F., 
Bertoldi, F., Chapman, N., Cieza, L., DeVries, 
C. H., Ridge, N. A., & Wahhaj, Z. 2004, ApJS, 
154, 396 



This 2-coluinn preprint was prepared with the AAS lATJiiX 
macros v5.2. 



19 



Table 1 
Calibration Coefficients and Uncertainties for BLAST05 



Band 



calib. coefF. 

[Xl0l2 Jy V-1] 



uncertainty Pearson correlation matrix 
[%] 250 Aim 350 Mm 500 ^m 



250 


7.61 


350 


3.16 


500 


1.56 



12 
10 



0.97 
1 



0.87 

0.96 

1 



Note. — Calibration coefRcients, calibration uncertainties, and Pearson 
correlation matrix, showing the relationship between errors in different 
bands for BLAST05. 



Table 2 
Flux Densities of BLAST05 Targeted Sources 



Name 


RA" 
[J2000] 


DEC" 
[J2000] 


Fli 
250^1111 


jx Density [Jy] 
350 Aim 


500 /im 


Pallas'" 








11.6 ± 1.5 


6.3 + 0.7 


3.7 + 0.3 


CRL 2688 


21'' 02'" 


'18! 75 


+36°4l'37''8 


113 + 14 


49 + 5 


21.6 + 1.7 


LDN 1014 


2l''24^ 


■"06= 


+49°59.'l 


22.4 + 2.7 


16.9+ 1.7 


8.9 + 0.7 


IRAS 20126+4104 
IRAS 21078+5211 
IRAS 21307+5049 


20'' 14" 
2l''09'' 
2l''32'' 


"25=1 
"25! 2 
■'31!5 


+41°13'32" 
+52°23'44" 
+51°02'22" 


590 ± 71 

1300 ± 160 

102 ± 12 


256 + 26 
590 + 60 

46 + 5 


97 + 8 
260 + 21 
18 + 1.5 


IRAS 22134+5834 
IRAS 23011+6126 
K3-50A'= 
K3-50C'= 


22''15'' 
23''03'' 
20''01'' 
20''01'' 


"09! 1 
"13!9 
"45!6 
"54! 2 


+58°49'09" 
+61°42'2l" 
+33°32'42" 
+33°34'15" 


418 ± 50 

95 + 12 

2100 ± 270 

1870 ± 240 


208 + 21 

59 + 6 

590 + 67 

1000 + 110 


82 + 7 
31 + 2.5 
270 + 26 
370 + 35 


W 75N'' 


20''38" 


"36! 5 


+42°37'35" 


4500 ± 540 


2000 + 200 


730 + 60 


DR21 {OH)"^ 
DR21'' 


20''39'' 
20''39" 


"00! 9 
"01!1 


+42°22'38" 
+42°19'43" 


9100 + 1400 
8800 + 1400 


4600 + 660 
3900 + 550 


2540 + 330 
920 + 120 


Mrk 231 


12'' 56"' 


'14! 23 


+56° 52' 25" 2 


6.0 + 0.8 


1.5 + 0.3 


... " 


Arp 220* 


15'' 34"' 


'57!21 


+23°30'09''5 


24.2 


9.8 


3.9 



Note. — Flux densities and associated uncertainties for BLAST05 targeted sources in the 250, 
350, and 500 fim bands. Quoted flux densities have been color-corrected. The uncertainties include 

the estimated contributions from calibration uncertainty as well as instrumental noise. 

^Positions arc nominal as given by SIMBAD. 

Pallas flux densities are average values from 4 observations taken at the following Julian Dates: 
2453534.26; 2453534.29; 2453535.22; and 2453535.40. 

■^Tlie K3-50 region is resolved into two sources by BLAST05. 

The W 75N region is resolved into three sources by BLAST05. 

''Mrk 231 was not detected at 500 /^m by BLAST05. 

Arp 220 is the absolute flux calibrator for BLAST05. Flux densities presented here should thus 
be considered as predictions based on a model fit to data from other instruments. Uncertainties in 
this model are given in § 14.11 
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Table 3 
Single Temperature SED Best Fit Parameters of BLAST05 Targeted Sources 



Name 


T 


/3 


5fir 


Distance 


^FIR 


Mdu.t 




[K] 




[W m'^] 


[kpc] 


[Lq] 


[Mq] 


Pallas" 


50.2 


0.2 ±0.2 


1.2 X 10"^^ 








CRL 2688" 


210 


0.4 ±0.2 


1.2 X 10"° 


1.25 


1.7 X 10^ 


3.2 X 10"^ 


LDN 1014 


12 ±3 


1.8 ±0.5 


(3.3 ± 1.2) X 10"^^ 


0.2 


0.38 ±0.14 


(3.7 ±3.5) X 10"^ 


IRAS 20126+4104" 


41 


1.1 ±0.2 


9.3 X 10"^^ 


1.7 


4.0 X lO'' 


0.67 


IRAS 21078+5211 


33.7 ±2.0 


1.0 ±0.2 


(1.2 ±0.1) X 10"^° 


1.65 


(1.0 ±0.1) X lO'' 


0.29 ±0.05 


IRAS 21307+5049 


32.9 ±2.5 


1.2 ±0.2 


(1.0 ±0.1) X 10"^^ 


3.6 


(4.3 ±0.5) X 10^ 


0.11 ±0.03 


IRAS 22134+5834 


23.5 ± 1.0 


2.1 ±0.4 


(2.0 ±0.8) X 10"^^ 


0.90 


520 ± 200 


0.11 ±0.07 


IRAS 23011+6126 


27.8 ±2.8 


0.6 ±0.2 


(4.3 ±0.5) X 10"^^ 


0.73 


71 ±9 


0.01 ±0.007 


K3-50A 


39.7 ±2.8 


1.7±0.2 


(6.2 ±0.5) X 10"^° 


8.5 


(1.4 ±0.1) X 10*^ 


8.1 ±1.4 


K3-50C 


9.3 ± 1.6 


3.9 ±0.5 


(2.3 ±0.8) X 10"^^ 


8.5 


(5.2 ±2.0) X lO'' 


(1.2 ± 1.1) X 10^ 


W 75N 


65" 


0.8 ±0.2 


(2.5 ± 1.3) X 10"" 


3 


(7.1 ±3.6) X 10^ 


1.0 ±0.1 


DR21 (OH) 


40"^ 


0.5 ±0.3 


(7.3 ±3.3) X 10"^° 


3 


(2.1 ±0.9) X 10^ 


4.9 ±0.7 


DR21 


29.2 ±3.0 


2.3 ±0.2 


(1.7 ±0.3) X 10"" 


3 


(4.9 ±0.8) X 10^ 


8.7 ±3.0 


Mrk 231 


43.9 ±2.5 


1.5 ±0.1 


(1.8 ±0.2) X 10"!^ 


175 000 


(1.6 ±0.1) X 10^^ 


(6.4 ±0.8) X lO" 


Arp 220"^ 


41.7 ±3.5 


1.3 ±0.1 


(5.9 ±0.8) X 10"^^ 


75 000 


(1.04 ±0.14) X 10^^ 


(7.8 ± 1.9) X lO** 



Note. — Parameters for a single temperature modified blackbody fit to the SED. T and /3 are the best fit to the BLAST05 and 
other data as indicated in the text. S-piu, is the total flux from the modified blackbody flt. Errors are from 100 Monte Carlo 
simulations of the flt. Distances given are from references indicated in § [s] FIR luminosity is base d on these adopted distances. 
MduBt is based on the formula A/dust — {S^D )/{k.Bi,{T)), with an assumed value of 10 for k fsee lChapin et al.||2008|) . 

^The data included in the flt do not accurately constrain the temperature, so no error bars are given. 

Data from [Harvey et al.| jlQTTD suggest a temperature of 65 K, which is used as a constraint in the flt. 

^Data from [Harvey et al.| ||1977[) suggest a temperature of 40 K, which is used as a constraint in the flt. 

Arp 200 is the primary calibrator for BLAST05. This is our best-fit to the other data, which we use to calibrate the BLAST 
photometry. 
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